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ABSTRACT 

In the article we present a hardware meteor detector. The detection principle 
is based on the electromagnetic wave reflection from the ionized meteor trail in 
the atmosphere. The detector uses the ANADIGM held programmable analogue 
array (FPAA), which is an attractive alternative for a typically used detecting 
equipment - a PC computer with dedicated software. We implement an analog 
signal path using most of available FPAA resources to obtain precise audio signal 
detection. Our new detector was verified in collaboration with the Polish Fire¬ 
ball Network - the organization which monitors meteor activity in Poland. When 
compared with currently used signal processing PC software employing real ra¬ 
dio meteor scatter signals, our low-cost detector proved to be more precise and 
reliable. Due to its cost and efficiency superiority over the current solution, the 
presented module is going to be implemented in the planned distributed detectors 
system. 

Subject headings: meteor engineering: general - meteor engineering: passive radio 


detection 
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1. Introduction 


One of the methods used for meteor detection is the radio meteor observation 


(McKinle y 1961 


Weitzen fe Ralston 1988 


Wislez 1995). The principle is based on the 


electromagnetic waves reflection from the meteor’s trail in the atmosphere, which is called 
the meteor scattering. When a meteor enters the atmosphere, it ionizes the electrons in 
atmospheric E-Layer (90-150km, Ionosphere layer) and leaves a plasma trail which, for a 
short time, works like a mirror for radio waves. Although the effect of the scattering is the 
strongest for 20-50 Mhz frequencies, it is possible to receive reflected signals from distant 
media radio stations operating in higher frequency range (90Mhz and above). This makes 
the phenomenon observable even by amateurs with a relatively low-cost equipment. The 
phenomenon of scattering is utilized extensively for meteor burst communication, which 
is a radio propagation mode allowing to establish short communication windows between 
distant stations. While the service of local media stations (like TV or FM) extends only 
slightly the optical horizon, reaching tens of kilometers, the maximum range the meteor 
scattering allows reception of signals from stations up to 2000km away ( Sugar 1964 ). The 
phenomenon is visualized in Fig. CD 


The number of meteoroids entering the atmosphere per hour changes periodically 
during the day. While in the morning we can observe the meteors approaching the Earth 
and the ones which are slower than our planet, in the evening only the meteoroids faster 
than the Earth enter the atmosphere. This fact causes characteristic diurnal variations of 
the meteor rate. The principle of this phenomenon is depicted in Fig. [2J 
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Fig. 1.— The principle of meteor scatter phenomenon. 
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Fig. 2.— The visualization of the Earth in meteoroids environment. 


In contrast to visual techniques, like video registration, which can not be carried 
out during the day, full moon or cloudy weather, the radio observations of meteors 
are of great interest as they can be performed continuously. The obtained count 
histogram is scientifically valuable and provides a number of meteor characteristics 
which require the analysis of the underdense and the overdense meteor distributions 


(jTrigo-Rodriguez et al 20081 ). We would like to note that the overdense and the underdense 
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meteors are the results of respectively heavy and light particles entering the atmosphere. 
The main coefficient - the mass index - provides information about the composition of 
meteoroids and their origins. While the high index is observed for showers from remnants 
of comets, the low index suggests the planetoid origins. The brief review of tasks, which 
may be accomplished employing the count histogram, is given below. 

1. The time occurrences of maximum of meteor showers are compared with the 
predictions to obtain the time-shifts whose origins have to be investigated. 

2. New meteor showers, observed in the form of abnormal activity, may be discovered. 

3. The overdense meteor mass distribution may be obtained by plotting logarithm T (the 
duration of phenomena) against the logarithm of the number of meteors with duration 
exceeding T. Only the phenomena longer than 0.5 seconds have to be included here. 
The slope should be 1 — s 0 , where s Q is the mass index of overdense meteors (see the 
left plot in Fig [3]). 

4. The mean mass index, which includes information from small and large particles, may 
be derived from the dependency between the logarithm of phenomenon amplitude 
and the number of detections above the given amplitude. The method of calculation 
of mean mass index s m is explained graphically on the right plot of Fig. |3j One can 
obtain required data in two ways: directly, by investigating detected signal amplitude, 
or indirectly, by utilizing parallel detectors set to various sensitivities. 


It has to be mentioned that by using a more sophisticated transmitter - receiver pair 
(active radar) some further properties of the meteoroids can be obtained. The array of 


receivers and high-speed signal corre’ 


ation computations can provide information about 


meteoroid velocity vector nKero 2012 ). Moreover, some observations of the ablation and 
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the fr agme ntation of the_ meteoroids are possible with a specially designed radar system 
( Malhotra & Mathews 2011). 



Fig. 3.— Mass index calculations. On the left - the overdense mass index estimation uti¬ 
lizing phenomena duration, on the right - the mean mass index estimation from amplitude 
histogram. 


In our project, we concentrate on a passive radio detection circuit for Polish Fireball 
Network (PFN) research activities. The PFN is t he organization which continuously 
monitors meteor showers and fireballs over Poland ([Wisniewski et ah 2012). The network 
consists of several base stations located all over the country, equipped with high-speed 
camera systems, recording and analyzing images in real time. They also use passive radio 
etection methods to enab l e observations du ring bad weather conditions or during the day 


( Zloczewski fe Szewc 2004 


Fictkiewicz 2012). In Fig. [4] we show an exemplary monthly 


radio count histogram provided by PFN, which includes Quadrantids meteor shower. For 
this purpose they used a PC computer, specially developed software (MetAN) and a 
Realistic DX-394 receiver tuned to distant, east European radio stations. 
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Fig. 4.— Monthly detection histogram from the Polish Fireball Network (January 2012). 
Light gray line - detec tions longer than 1 s, black line - detections longer than 0 s (all 


detections) flPFN 20131 1 .The peak refers to the Quadrantids meteor shower on 3 rd January 


2012. The measurements were binned in 0.1 solar longitude time intervals (about 2.43 hour). 


The PFN PC computers running 24 hours a day and detecting meteor radio signals are 
neither energy efficient nor portable solutions. It is to be noted that while a single station 
is still able to provide useful information about meteors, a system of multiple detectors, 
located at different sites and directed to various distant stations, may significantly improve 
the measurement efficiency and extend the analysis capabilities. These two main issues are 
discussed in detail below. 

1. Firstly, we should note that the detector robustness depends on local electromagnetic 
environment. The sites far from strong radio stations are generally preferred since 
there are fewer false detections resulting from interfering frequencies. However, the 
best sites usually suffer from lack of access to the power, thus ideally the detection 








module should be self-sufficient e.g. powered by solar cells. 


2. The reception of reflected signal depends on the relative position of the receiver, 
meteor trail and remote radio station. By directing the detector toward different 
stations we would gain additional information about the brief position of trail 
relatively to the detector. Moreover, a distributed system of detectors, spread over 
hundreds of kilometers, should provide far better meteor localization capabilities and 
would work like a radar. We plan to investigate this issue in the future. 


In the article we introduce a dedicated detection electronic circuit, easy connectible 
to any radio receiver, which can work as a module in a remote, battery powered, meteor 
scatter detector. For our task we selected a field programmable analog array (FPAA) 
manufactured by ANADIGM. We also intend to show that for such a purpose an analogue 
array is an interesting alternative to commonly used digital signal processing units. 


The circuit is dedicated to planned distributed system of remote meteor detectors. 
We also consider the amateur application of our module. With a growing number of 
astro-amateurs and dedicated self-sufficient equipment, it may be possible to explore the 
meteor phenomena as it is performed. for variable stars in the American Association of 


Variable Star Observers (AAVSO) ( Williams fe Saladvga. 2011 ). 


The article is organized as follows. In Chapter 2 we introduce the internal structure of 
the ANADIGM FPAA circuit and its principle of operation. Our signal chain implemented 
for meteor detection is explained in detail in Chapter 3. Finally, in Chapter 4, we describe 
our verification methodology based on the comparison with the results obtained by MetAn 
software, which is currently used for meteor scatter detection by the Polish Fireball 
Network. The results of the tests are presented and analyzed in Chapter5. We summarize 
the results of our research in Chapter 6. 
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2. ANADIGM Field Programmable Analog Array 


The AN231E04 device represents the th ird generati on of the dynamically reconfigurable 


FPAA AnadigmApex offered by Anadigm fires, doc, ll ). It is powered by 3.3 V supply 


voltage. The previous FPAA family, called Anadigm Vortex, was powered by 5V, but its 
DC performance was much worse. 


The FPAA circuit consists of four configurable analog blocks (CAB) built on the 
switched-capacitance principle, surrounded by a network of programmable connections. 
There are also seven programmable input/output cells on the silicon chip. Three of 
them can have access to a specialized chopper amplifier resource, which allows accurate 
amplification of very low energy input signals. All internal signal paths of the block are 
differential. 


10 cells contain both passive and active circuitry, which allows direct signal input and 
output, building active filters, sample and hold circuits, digital inputs, and digital outputs. 
The frequency response of continuous time input and output filters is determined by a 
combination of internal programming and external components. 

The block diagram of the AN231E04 device is shown in Fig. [5l where one can see its 
basic elements: configurable analog blocks (CAB), input and output cells, look up table 
(LUT), clocking circuits, communication interfaces and reference voltage sources. 

The new AN231E04 devices can be reconfigured dynamically (on the run), which 
enables e.g. auto ranging, auto calibration and multimode work (allows to build many 
functions in a single IC). The microcontroller connected to the AN231E04 device can load 
the new device configuration data to the device still working in its previous configuration. 
The configuration data can then be swept in a single clock cycle. This new and prospective 
feature of the dynamic reconfiguration ability of the AN231E04 devices enables the creation 
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of a modern analog system, which can be (completely or in part) rebuilt in real time. 



Fig. 5.— Internal structure of ANADIGM field programmable analogue array. 


We chose this particular FPAA because of its proved very good cap ability of audio 


signal processing (Falkowski & Matcher 2 011 


Falkowsk i & M alche r 2010). By using 


different configurations and changing individual parameters we were able to create a circuit 
well suited for our purpose. Such designing flexibility makes FPAA a very attractive 


prototyping device. Moreover, the FPAA devices ha ve bee n recent 


detectors of QRS complexes in the EGG signal (jMalcher et al. 2010 


y used as h ardware 


Malcher et al. 2009), 


which is an application similar to our audio meteor detector. 


In comparison with popular Field Programmable Gate Arrays (FPGA), which also 
might be utilized as a meteor detector, FPAAs have several advantages. Below, we 
enumerate the most crucial ones for our project. 


1 . 


FPAA devices introduce much lower electromagnetic disturb ance, w hich should be 
considered when analyzing weak analogue signals (j Mocha fc Woznica 20121 ). 


2. The implementation process is very simple and efficient as the construction of 
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signal path may be accomplished with an easy-to-use, drag-and-drop, diagram based 
environment - AnadigmDesigner (FPGAs have to be programmed using dedicated 
VHDL or Verilog register transfer languages). 

3. Nearly all modifications of the existing signal chain may be tested by changing single 
values in user-friendly environment (the properties are clear and simple, like filter 
bandwidth or gain). 

4. The analogue filters implemented in FPAA are more reliable and efficient in contrast 
to their digital versions (the infinite response HR digital filters suffer from low 
stability while the finite response FIR implementations, due to the high number of 
parallel calculations, usually require specialized high-efficiency units, like digital signal 
processors DSP ). 

5. The FPAAs do not require any analogue-to-digital conversions since they operate 
in the analogue domain, therefore the number of elements and the total cost of the 
circuit is lower. 


3. Detection circuit 

Our signal chain consists of several blocks (see Figj6]): primary band-pass filter 10Hz 
- 600Hz (BPF), full-wave rectifier (RECT), secondary filters (LPF1, LPF2+RCLPF), 
primary comparator (COMP1) and analogue realization of the time off-delay function 
(TOFF). We used nearly all available analogue cells of ANADIGM FPAA. 
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Fig. 6.— FPAA implementation of the meteor scatter detector: INPUT - input audio signal, 
BPF - 10-600Hz band-pass filter, RECT full-wave rectifier with a low-pass filter (envelope 
detector), LPF1 - low-pass filter 1 (20Hz), LPF2+RCLPF - low-pass filter 2 (0.1Hz), TR 
- threshold level adjustment, COMP1 - comparator 1, PEAK DET - peak detector, REF - 
detection dead time reference level adjustment, COMP2 - comparator 2, RAW OUT - output 
of the detector without time off-delay circuit, TOFF OUT - output of the detector with time 
off-delay circuit. 

As the analyzed signal is in the audio frequency range, we decided to choose a primary 
band-pass filter 10 - 600 Hz. It is due to the fact that most of the spectral components are 
within this range. As a proof, we present a plot of a part of registered radio signals in the 
time-frequency domain in Fig. [71 The corresponding audio signals are the representative 
examples obtained by the PFN radio meteor station in our experiment. We depicted various 
phenomena, with long and short duration, high and low amplitudes. We checked that all 
investigated audio signals in our recording were well within the assumed frequency range. 
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Fig. 7.— Audio signals of four meteor scatter phenomena in the time-frequency domain. The 
examples include various amplitudes and different durations. By dashed line we indicate the 
cut-off frequency (600Hz) of our detector. 

The secondary filters (LPF2 with RC LPF) provide information about the mean level 
of the signal power. We constantly compare the output of 20 Hz filter LPF1 with the biased 
(TR - threshold) output of the 0.1 Hz filter LPF2. The biasing prevents the false detections 
induced by temporal noise fluctuations. Both filters outputs represent the level of signal 
power: the first one detects rapid level changes, the second one shows a long-term trend. 
The first comparator (COMP1) is set high when the LPF1 output signal exceeds the LPF2 
output. COMP1 output is our first-stage detector (RAW OUT). 

However, we observed that during a single radio station signal appearance the RAW 
OUT output had several positive detections, although there should be only one meteor 
phenomenon detected. A solution to this problem is provided by an analogue version of a 
time off-delay (TOFF) circuit. It is realized with a filtered peak detector (PEAK DET), a 
reference signal and a comparator. When the positive edge of the first comparator (COMP1) 
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appears, the output of the peak detector (PEAK DET) and the output of the whole device 
(TOFF OUT) immediately rise and are held as long as the output of the first comparator 
(COMP1) remains high. When it switches to the low state, the output of the peak detector 
(PEAK DET) slowly falls down as an RC circuit in the peak detector discharges. As long as 
the reference (REF) is not reached, the second comparator (COMP2) remains high. When 
the reference level is reached, the second comparator (COMP2) changes the output state 
to low, but if there is another detection in the meantime, the device output is not affected 
and remains high. Adjusting the reference (REF) or the RC time constant in the peak 
detector can change the time delay. We decided to choose a 3 second interval, which is also 
our estimation of the minimal interval between the subsequent meteors. The disadvantage 
of this method - the time delay of the detection end - can be easily overcome by reducing 
measured length of the phenomena by a corresponding and constant time interval. We 
would like to note that the detector will produce an additional detection if there is a longer 
than 3-second silence in the radio program. It may be easily mitigated by extending the 
delay; however, it also reduces the detector ability to distinguish close events. 

The principle of operation is presented in Fig. |8j We used two meteor scatter audio 
signals - a long one (on the left) and a short one (on the right) - and showed the most 
important waveforms obtained inside the designed FPAA circuit. The figures present the 
screenshots registered by the TECTRONIX TDS2014B digital oscilloscope. The detections 
start when the signal from LPF1 output exceeds LPF2 output. The operation of the 
analogue time off-delay circuit is clearly visible on the PEAK DET OUT and the TOFF 
OUT waveforms. 
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Fig. 8.— An example of signals generated inside FPAA during a very long (over 15 seconds, 
left figure) and a short (several hundreds milliseconds, right figure) meteor scatter signal. 


4. Verification methodology 

We used two-hour recording provided by the PFN. Their equipment for our tests 
included a 5-element log-periodic antenna (Fig. [9]), 9dBi gain, provided by VPA-SYSTEM. 
The antenna was directed toward a distant Belorussian transmitter. We provide all 
necessary transmitter’s information in Tab. [0 An audio signal from the Realistic Pro and 
DX-394 receiver tuned to 70.10 MHz was captured and saved on a PC computer as a wav 
audio hie. Since the Polish broadcast FM radio stations use the frequency band of 88-108 
MHz, none of them should affect the measuring meteor scatter audio signals and cause false 
detections. 
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Tabic 1: Transmitter characteristics. 


location 

Mahilyow/RTPS Polykovici 

longitude 

30° 19’ 36.41” 

latitude 

53° 59’ 25.92” 

frequency 

70.10 Mhz 

program 

BR Radio Mahilyow 

power 

4.00 kW 

azimuth from receiver 

67° 

distance from receiver 

667 km 



Fig. 9.— 5-element log-periodic antenna used by the Polish Fireball Network (direction: 
north-east, azimuth 67°, horizontal polarization, height 3m). The picture was taken from 
below the antenna mast. 

As the detection algorithms implemented in the MetAn software and in our detector 
are highly dependent on a threshold level, we decided to repeat the detection procedure for 
several values of the threshold. We did it for both the MetAn software and our hardware 
detector. For this purpose we used several parts of previously mentioned recording, to 
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reduce the measurement time for each sensitivity level. The constructed, well defined, 
test signal included 40 clearly audible signals extracted from 2-hour recording. They were 
combined with the proceeding and following 10 seconds, thus our test recording was about 
15 minutes long (instead of 2 hours). 

The resulting dependencies of the number of detections and the total detections 
duration are presented in Fig. [TQ1 It is visible that each resulting curve shows a 
characteristic break related to the rising number of false detections due to the noise level 
comparable to the threshold. We decided to set both detectors to their maximum and 
medium sensitivities, which are represented by the vertical solid and dashed lines in both 
corresponding figures. 

Firstly, we selected the maximum sensitivity thresholds, for which the number of 
obtained detections was as close as possible to the number of given phenomena (40). It 
also let us estimate the thresholds that should not be exceeded to avoid a rapidly rising 
number of false detections. The medium sensitivity thresholds were chosen for the number 
of detections equal to about 30. It corresponds to the situation, where the detectors have a 
slightly higher detection threshold and consequently they miss about 10 weakest signals. 

We would like to note, that the total detections duration for MetAn software is 
significantly larger because we did not know the additional time interval added after each 
detection. In the case of our hardware detector the delay of detection end was known - 
exactly 3 seconds - thus we were able to remove this additional offset. 

Finally, the main recording (2 hours long) was given at the FPAA input and at the 
line-in input of a PC computer which was running the MetAn software detector. We 
observed simultaneously the output of both detectors and the input signal. The detections 
for both medium and maximum sensitivities were compared. All events detected by both 
detectors were carefully examined using headphones to check if the audio signal was audible. 


18 


In this way we classified detections as the false and the true ones. 



Fig. 10.— Number of detections (dashed black line) and total detections duration (solid 
gray line) as a function of the threshold level for the MetAn software (on the left) and for 
our FPAA based detector (on the right). By the vertical solid line we suggest the threshold 
level for maximum sensitivity. By the vertical dashed line we suggest the threshold level for 
medium sensitivity. 


5. Results 

The results are presented in Fig. [HI Each event time can be computed by adding 
TimeA and TirneB. The total duration of our recording (2 hours) was reduced to 6500 
seconds in Fig. [Has the last part did not include any detection. We chose full circles and 
squares to identify false detections of respectively FPAA and MetAN detector, empty circle 
markers for positive FPAA detections and crosses for MetAN positive detections. 


















19 


6000 
5500 
5000 
4500 
4000 

Tv 3500 

CD 

qj 3000 

E 

P 2500 
2000 
1500 
1000 
500 
0 

Tima A [s] 

6000 

5500 

5000 

4500 

4000 

TT 3500 

CO 

0) 3000 

E 

P 2500 
2000 
1500 
1000 
500 
0 

10 

Fig. 11.— Time plot of detections for medium (upper plot) and maximum (lower plot) 
sensitivities. The event time can be calculated by adding TimeA and TirneB. 
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The results for both sensitivities show that many of the detections overlap. Almost 
all detections of the MetAN software were also detected by FPAA (except for 4 and 2 
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detections respectively for maximum and medium sensitivities), which confirms that FPAA 
covers the domain of MetAN. However, we would like to indicate that our detector found 
a lot of additional events in the full recording, which were not recognized by the software 
detector (red circles without blue crosses inside, respectively 22 for maximum and 8 for 
medium sensitivity). Moreover, the number of false positive detections is higher for the 
MetAN analyzer, thus our detector is more resistant to noise fluctuations, which may 
influence the detector’s decision. We summarize the results in Table [2] 

The quality of measurement, hence the quality of obtained results, depends on both the 
number of true detections (the signal) and the number of false detections (the noise floor). 
Ideally, the detector should maximize the signal level while keeping the low noise floor (i.e. 
maximize the signal-to-noise ratio). The results of our experiment for both sensitivities 
prove that this done much better in our solution than in the competitive detector. 

Table 2: Detections summary for both detectors and for both sensitivities. 



Medium sensitivity 

FPAA MetAN Common 

Maximum sensitivity 

FPAA MetAN Common 

True positive 

29 

23 

21 

55 

37 

33 

False positive 

0 

0 

0 

1 

3 

0 


6. Conclusions 

In the article we presented a hardware meteor radio detector which can be attached 
to any radio device to detect the meteor scatter phenomena. The signal processing 
circuit of our detector was implemented in ANADIGM field programmable analog array 
(FPAA) using nearly all its available resources. A signal path contains a rectifier, a peak 
detector, several comparators, low-pass and band-pass filters. The output signals of the 
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filter provide information about signal power in the audio frequency range. Additionally, 
the implemented time off-delay analogue circuit prevents multiple detections during long 
scatter phenomena. Such an FPAA device enabled flexible signal flow designing to fulfill 
the detector requirements. 

We compared the efficiency of detection of our implementation with the current 
solution used by the Polish Fireball Network (PFN) - the MetAN software dedicated to a 
PC computer. We presented the results of the analysis of 2-hour recording received from 
a PFN meteor station receiving Belarusian radio signals. By setting both detectors to 
medium and maximum sensitivities we compared the effectiveness of both solutions for 
different detection threshold levels. 

The results proved a significant difference between the detectors. The FPAA detector 
proved its ability to pick up noticeably more audio signals for both thresholds. When set to 
the maximum sensitivity, it detected even hardly audible scatter events while maintaining 
its low false detection rate (only 1 false positive per 57 detections). In comparison, current 
software detector found significantly fewer events and had lower true detection percentage. 
Our experiment confirmed that the measurements employing FPAA detector are of higher 
efficiency, thus the outcomes have a better signal to noise ratio. 

Currently, we are developing a ready-to-go electronic board with peripheries 
incorporating the presented FPAA-based detector, recording on an SD-card and an optional 
screen presenting counts histograms. Such a device could be easily handled not only by 
professional meteor researchers, but also by amateur astronomers who wish to contribute 
to additional radio observations. With its efficiency and very little power consumption 
(comparing to a PC computer), our detector might replace current software-based detecting 
methods. We also wish to set a network of remote detectors which will provide complex 
data and enable further investigation of its radiolocation capabilities. 
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